Random triangle in square: geometrical approach 
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The classical problem of mean area of random triangle within the square 
is solved by a simple and explicit method. Some other related problems are 
also solved using Mathematica. 

I. INTRO 

We call our approach geometrical as instead of considering 6-fold integral in abstract 
space we consider random triangle (RT) inside the plane rectangle when all possible cases 
are explicitly apparent. 

Area of triangle with vertices pl=(xl,yl), p2=(x2,y2), p3=(x3,y3) is equal to 
s = §(xl(y2 - y3) + x2 (-yl + y3) + x3 (yl - y2)). (1) 

Let points pi, p2, p3 are randomly (with constant differential probability function) dis- 
tributed over the rectangle with sides A, B. What is the mean area of triangles with vertices 
pl,p2,p3? 

Answer is evident: zero, as any given triangle corresponds to 6 cases of full permutation of 
three points at vertices of the triangle. Mean area of this 6 triangles, as given by (1), is zero. 
But if we take triangle as geometrical figure and if we consider an area of such a figure as 
positive value, then we must take absolute value of s in formula (1) and... calculation of 
relevant integrals become impossible even for Mathematica. So Michael Trott in his recent 
brilliant paper in Mathematica Journal [1] found 496 different integrals each over subregion 
with the same sign of s, and then used Mathematica to solve such an enormously difficult 
task. Needless to say that M.Trott's stunning skill of using Mathematica is far out of scope 
of ordinary reader (as me, e.g.), so I've spend some three weeks in searching a more simple 
solution. The result is most easily get by the explicit geometrical approach. 

II. GEOMETRICAL APPROACH 

Here we consider RT in square (=right rectangle) with side length A. First observation 
is that due to points symmetry we may simplify problem by considering particular relation 
between points. For example, as we will do here, we may consider only case xl < x2 < x3 
with due account of normalizing condition. 

First point (pi) may take any position inside the square, so the region of integration over 
xl, yl is < xl < A; < yl < A at all cases considered further. 

Now we should discriminate two cases of relation between ordinates of 1st and 2nd points: 
1) y2 > yl and 2) y2 < yl. 
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A. y2 > yl 



In this case, important is relation between two angular coefficients kl and k2: 
kl — (A — yl)/{A - xl); kl = (y2 - yl)/(x2 - xl). 

k2 > kl. In this case the line (pi, p2) crosses the upper side of square at the point (x31m, A) 
with x31m = (A - y2)/k2 + x2, see Fig.l , panel 1A. 
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FIG. 1. Ten different cases of integration regions, bounded by solid lines 

The region of integration over x2, y2 is: 
xl < x2 < A; y2\m < y2 < A, y2\m = kl(x2 - xl) + yl. 

The first region of integration over x3, y3 is : x2 < x3 < x31m; k2{x3 — x2) + y2 < y3 < A, 
region 1+, panel 1A, Fig. I. 

At this region formula (1) gives positive value as points (pl,p2,p3) make right-handed 
system: moving in direction (pi — > p2 — > p3 — > pi) we make counter-clockwise rotation. 
Now we are ready to calculate first integral: 

r-A r-A rA r-A r-xZlm r-A J\£> 

II = dxl dyl / dx2 / dy2 / dx3 / (s)dy3 = —r^- 

JO JO Jxl Jy21m Jx2 J k2(x3-x2)+y2 34560 

Second integral appears from region "under" the first integral's region, region 2-, panel 1A, 
Fig. 1, and here formula (1) should be taken with negative sign: 

fA rA rA rA rx3lm rk2(x3~x2)+y2 23A S 

12 = dxl dyl / dx2 / dy2 / dx3 / (~s)dy3 = —r—. 

Jo Jo Jxl Jy21m Jx2 Jo 34560 

Note that 2nd integral differs from 1st one only by integration boundaries over y3 (and by 
sign of s). 
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Also, the interesting exact relation occurs between numerical values of two considered inte- 
grals : 12 = 23/1. 

Last integral at the case k2 > kl, corresponding to region 3-, panel 1A, Fig. 1, is; 

[A [A [A [A rA [A 748 

/3 = / dxl dyl / dx2 / dy2 / dx3 / {-s)dy3 = — — = 140/1. 

JO JO Jxl Jyllm JxZlm Jo 1728 

Note that all three multiple integrals have the same first four particular integral regions and 
we may write down them in a more compact form, but we will not do this pure "decorative" 
operation. 

k2 < kl. In this case and still at yl < y2, the line (pl,p2) crosses the right side of square; 
now the region of integration over y2 is yl < y2 < y21m, and we have two different regions 
of integration over p3: 

x2 < x3 < A, k2(x3 — x2) +y2 < y3 < A, with positive s, region 4+, panel IB, Fig. l,which 
gives 14, and 

x2 < x3 < A, < y3 < k2(x3 — x2) + y2, with negative s, r. 5-. p. 1A, Fig. 1, which gives 
/o. 

Therefore we have two additional integrals: 



r-A r-A r-A ry21m r-A rA 194 8 

14 = dxl dyl dx2 dy2 dx3 (s)dy3 = — = 19/1. 

Jo Jo Jxl Jyl Jx2 Jk2(x3-x2)+y2 34560 

rA rA rA ry2\m rA rk2(x3-x2)+y2 37A 8 

15 = dxl dyl / dx2 / dy2 / dx3 / (s)dy3 = — = 37/1. 

Jo Jo Jxl Jyl Jx2 Jo 34560 



B. yl <yl 

Now important is relation between two coefficients k3 and k4 : 
k3 = yl/(A - xl); k4 = (yl - y2)/(x2 - xl). 

k4 < k3. If k4 < k3, then the line (pl,p2) crosses the right side of square , see panel 1C, Fig. 
1, and we have the case completely analogous to the case considered in the previous section 
(see also panel IB) and two integrals, 16 with positive s and 17 with negative s, are equal to 
14 and 15 respectively. Here our geometrical approach is particularly explicitly demonstrate 
his power: it is sufficient to look at panels 1A - ID of the Fig. 1 to be convinced that 
actually we have only two different cases, one case when line (pl,p2) crosses two opposite 
sides of square and second case when line (pl,p2) crosses two adjacent sides of square. 
k4 > k3. Now the line (pl,p2) crosses lower side of square, panel ID, Fig. 1, and we have 
in essence the case coinciding with previous one, panel 1A, Fig. 1, so we have another three 
integrals with values actually found before: 18=11, 19=12, and 110=13. 
Now sum of all 10 integrals is equal to // = llA 8 /864. The normalizing coefficient is found 
by calculating a sum of abovementioned 10 integrals with (+/-s) = l in all cases which gives 
J J = A 6 /6, that is 1/6 of volume of hypercube with side A. So the mean area of random 
triangle inside the square is II/ J J = 11/144 of host-figure's square, A 2 . 
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III. RT IN RECTANGLE 



Our geometrical approach allows easily to calculate also the mean area of random triangle 
when host-figure is rectangle. Being experienced with the square case we consider here only 
two cases leading to 5 integrals. 

As result, we present a simple and transparent Mathematica's code for calculating the mean 
area of random triangle in rectangle with sides A and B. 

(* y2\,>\,yl *) 

kl=(B-yl)/(A-xl) ;k2=(y2-yl)/(x2-xl) ; Y2=kl* (x2-xl)+yl ; Y3=k2* (x3-xl) +yl ; 

(* k2\,>\,kl *) 

X=(B-y2)/k2+x2; 

Il:=Integrate[s,xl,0,A,yl ,0,B,x2,xl ,A,y2,Y2,B, x3,x2,X,y3,Y3,B] ; 
(* Il=A~4*B~4/34560 *) 

12 : =Integrate [-s , xl , , A , y 1 , , B , x2 , xl , A , y2 , Y2 , B , x3 , x2 , X , y3 , , Y3] ; 
(* 12=23*11 *) 

13 : integrate [-s,xl,0, A, yl,0,B,x2,xl, A, y2,Y2,B, x3,X,A,y3,0,B] ; 
(* 13=140 II *) 
(* k2\,<\,kl *) 

14 : =Integrate [s , xl , , A , y 1 , , B , x2 , xl , A , y2 , y 1 , Y2 , x3,x2,A,y3,Y3,B] ; 
(* 14=19*11 *) 

I5:=Integrate[-s,xl,0,A,yl,0,B,x2,xl,A,y2,yl,Y2, x3,x2,A,y3,0,Y3] ; 
(* 15=37*11 *) 

115=11+12+13+14+15; (* I15=11*A~4*B~4/1728 *) 

(*calculation of normalizing coefficient*) 

s=l;Jl=Il;(* Jl=A~3*B~3/432 *) ; s=-l ; J2=I2 ; (* J2=5*J1 *) 

s=-l; J3=I3; (* J3=18*J1 *) ; s=l ; J4=I4; (* J4=5*J1 *) 

s=-l; J5=I5; (* J5=7*J1 *) ; J15=J1+J2+J3+J4+J5; (*J15=A~3*B~3/12*) ; 

MeanSquare0fRandomTriangleInRectangle=I15/J15; (* = (11/144) *A*B *) 

(*MeanSquareOf RandomTri angle I nRect angle = 11/144 of rectangle's square*) 



IV. RT IN SQUARE FRAME 

Now we consider the related problem of random triangle in square frame. Let three 
points are randomly (with constant differential probability function) distributed along the 
sides of unit square (side length and square being 1). What is the mean area of triangles 
formed by these points as vertices? 

The solution is elementary, but we consider it for pure pedagogical purposes. First obser- 
vation is that due to symmetry of square and due to symmetry of points it is sufficient to 
assume that 1st particle is at bottom side of square. Then four different cases should be 
considered. 

2nd particle is also at bottom side 

Let 3rd particle moves with constant linear velocity along all sides of the the square. We 
are looking for value of integral over coordinates of 3rd particle and then over coordinates 
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of 2nd particle which is allowed to move only along the bottom side of square. 
In the Mathematica's language we should calculate the following path integral: 



s[xl_,yl_,x2_,y2_,x3_,y3_] :=(l/2)*Abs [(xl*(y2-y3)+ 
x2*(-yl+y3)+x3*(yl-y2))] : 

II :=Integrate[ (Integrate [s [xl ,0,x2,0 ,x3, 0] ,{ x3,0,l}]+ 
Integrate [s [xl, 0, x2, 0, l,y3] ,{ x2, 0,1}, {y 3,0,1}]+ 
Integrate[s[xl,0,x2,0,x3,l] ,{ x2, 0,1}, { x3,0,l}] + 
Integrate [s [xl, 0, x2, 0, 0, y3] ,{ x2, 0,1}, { y3,0, 1}] ) ,{x2,0, 1}] 

Result is II = \ — xl + xl 2 . 

2nd particle is at the right side of the square 

The relevant integral which we do not write down is equal to 12 = u ~ Sx ^+ 3xl ■ 
2nd particle is at the upper side of the square. 13 = u ^ x ^+ 6xl2 , 
2nd particle is at the left side of the square. I A = 6 + 2 ^+ 3xl2 , 

Sum of these for integrals gives 114 = II + 12 + 13 + H = ^ - 2x1 + 2xl 2 . Now 

Integrate [1 14, {xl, 0, 1}] gives 5/2. The normalizing coefficient is evidently 16, as we 

calculate 4 path integrals (over 3rd particle) which of them has path length equal to 4. 

Final result is: the mean area of random triangle inscribed in unit square is 5/32. 

Let us look at this value (and check it!) from another point of view. We divide all sides 

of unit square to 10 equal parts and let each of three particles take all mid-points of these 

40 parts. Then we have 40x40x40x40=640,000 triangles with mean area (as calculated by 

Mathematica) equal to 249/1600=498/3200 that is very close to 5/32. 

So we understand more vividly in what sense the mean area of random triangles inscribed 

in the square is 5/32. 

It is very interesting to compare these two "mean" values 11/144 and 5/32. First value is 
22/45, that is almost exactly 1/2, of the second one. That is mean square of triangles with 
vertices randomly distributed all over the square is almost exactly half of mean square of 
triangles vertices of which are allowed to occur only at sides of square. 
The reason of considering this last problem is originally related to my attempts to find 
simple solution of first problem. Is seemed to me that by solving the problem of inscribed 
triangles I could somehow solve the first problem also. Some hazy ideas about differentia- 
tion/integration connection between two problems unfortunately gave no yield and I solved 
these two problems separately. 

What is left is the problem of mean volume of tetrahedron in the cube (M.Trott, personal 
communication ). I hope that geometrical approach will also help in this much more difficult 
problem. But if the geometrical approach managed to reduce the number of integrals from 
496 cumbersome ones in original solution by M. Trott to 5 very simple integrals, hopefully 
it will help in tetrahedron-in-cube problem as well. 

Numerical value get by Mathematica gives 1/72, but this is not exact value, this is value 
get by Mathematica's command Rationalize [NumericalValue, 10~-4] . 
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